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Abstract. We investigate the field- induced insulator-to-superfluid transition of 
bosonic quasiparticles in S = 1/2 weakly-coupled dimer antiferromagnets. In presence 
of realistic disorder due to site dilution of the magnetic lattice, we show that the 
system displays an extended Bose-glass phase characterized by the localization of the 
hard-core quasiparticles. 
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1. Introduction 

The quantum phase transition from the insulating to the superfluid state in strongly 
correlated bosonic systems is realized in a number of model systems, ranging from 
Josephson-j unction arrays [l] to cold bosons in optical lattices j2j. Recently it has been 
realized that quantum antiferromagnets also offer a clear example of such transition 
and this has stimulated an intense experimental and theoretical activity on the subject 
[H El E] ■ In the context of quantum magnets the bosonic degrees of freedom correspond 
to the local deviations from a given ground state, either ordered or disordered, and their 
condensation corresponds to the occurrence of spontaneous antiferromagnetic long-range 
order breaking a planar rotational symmetry. 

A particularly intriguing question in the context of bosonic quantum phase 
transitions is their fate in presence of lattice disorder. General considerations [Zj and 
specific examples coming from microscopic Hamiltonians |H] show that a novel Bose- 
glass phase is induced by disorder between the superfluid and the insulating phase, and 
it corresponds to a localized phase for the bosonic degrees of freedom. Nonetheless, 
the experimental observation of such a phase has been so far elusive. It appears 
therefore tempting to imagine its realization in quantum magnets, where disorder can 
be introduced in a highly controlled manner through site-dilution of the magnetic lattice 
|Hj and a rich variety of experimental probes are available to detect the specific features 
of this phase. 

In this paper we show the emergence of a Bose-glass phase in a site-diluted S — 1/2 
Heisenberg bilayer in a magnetic field. After a general description of the superfluid- 
insulator transition(s) in quantum magnets, we show how to recast the site-diluted 
bilayer in terms of a disordered Bose-Hubbard model for hard-core bosons. Then we 
present quantum Monte Carlo data on the original spin Hamiltonian, clearly supporting 
the scenario of an extended Bose-glass phase in a realistic quantum magnetic model. 

2. Superfiuid-insulator transitions in quantum magnets 

Although quantum magnets are nearly perfect Mott insulators in terms of electronic 
properties, they admit a well-known bosonic representation in terms of local elementary 
spin deviations with respect to a given spin configuration. In the case of classically 
ordered magnets, such spin deviations are the so-called magnons, namely bosonic 
quasiparticles carrying a spin 1. A perfectly ordered magnet, such as a Heisenberg 
ferromagnet in its ground state, represents therefore a trivial insulating state (the 
vacuum) for such quasiparticles. More complex ordered ground states, in which none 
of the spin components reach its saturation value, can be mapped onto bosonic states 
with a macroscopic number of bosons displaying long-range phase coherence, namely on 
superfluid bosonic states. The simplest example of such a state - and the one relevant for 
the remainder of the present paper - is that of a quantum Heisenberg antiferromagnet 
in a uniform magnetic field, displaying canted antiferromagnetic order. If the field is 
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applied along the z axis, using, e.g., the Villain spin-boson transformation ^U] the 
deviation of the z spin component from its saturation value at site i, S — S Z} i, represents 
the local number of bosons Hi, and the trasverse spin components S i ~ y/S exp(±j</>j) 
(in the limit of large spin S) bring the information on the local phase fa. Long- 
range antiferromagnetic order transverse to the field corresponds to long-range phase 
coherence of the bosons [TT]. The transition from a fully ferromagnetic state to a 
canted antiferromagnetic state is realized for instance in a spin-5 1 quantum Heisenberg 
antiferromagnet on a hypercube of dimensions D > 2 when an applied uniform field 
h = gUbH/ J (in units of the Heisenberg exchange coupling J) is driven below the 
critical value h c = ADS. This transition can be conveniently described as a band- 
insulator-to-superfluid transition in the bosonic language. 

A more complex bosonic insulating state is realized by quantum magnets having 
a quantum-disordered ground state. This is the case of systems of S = 1/2 
antiferromagnetic dimers with intradimer coupling J, coupled together through weaker 
antiferromagnetic couplings J' < J, with general Hamiltonian 

H = jJ2' s i-S j + J'^SrS m -hjY^St . (1) 

(ij) (lm) i 

where the (ij) nearest-neighbor bonds are intra-dimer bonds and the (lm) bonds are 
inter-dimer ones. In the limit of J' = each dimer is obviously in a singlet state 
\s) = I TO - I lt))/V2 with total spin Stot — 0. This state naturally represents the 
vacuum for three S'tot = 1 bosonic triplets \t±) = | ft), | II), and \t ) = | TD + I It))/v2, 
separated by a gap of J from the singlet ground state. A weak coupling J' between the 
dimers slightly perturbs the simple picture of dimer singlets, and the state of each dimer 
is properly described as an incoherent mixture of the singlet state with the triplet states, 
retaining the full SU(2) symmetry of the Hamiltonian: the pure-state probabilities, 
corresponding to the diagonal elements of the dimer reduced density matrix 
p(\ s )) ^ P(l^o)) — p(\t+)) — p(\t-))i an d all off-diagonal terms are vanishing. This 
picture of a so-called dimer-singlet ground state is valid below a critical value of the 
ratio g = J'/ J, strongly dependent on the dimensionality of the resulting lattice and on 
the specific geometry of the weak couplings. In what follows we assume the system to 
be well inside the dimer-singlet regime in zero field. 

Applying a magnetic field to each isolated dimer (J' = 0) along the quantization 
axis, the two triplets \t~) and \t°) remain well separated from the singlet ground state, 
while the triplet is brought to degeneracy with the singlet for a field h — 1, at 
which the system jumps to a fully magnetized state. In presence of a finite interdimer 
coupling J' > the application of a weak field h < 1 reduces the symmetry of the 
Hamiltonian from SU(2) to U(l), and correspondingly in the dimer-singlet phase the 
reduced density matrix p( dim ) for each dimer has p(\to)) ^ p(\t + )) = but the last 

equality guarantees that the system does not develop a finite magnetization along the 
field, so that the ground state has a higher symmetry than the Hamiltonian. Moreover 
the mixture remains incoherent to mantain the U(l) symmetry. 
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Nonetheless, the interdimer J' interaction is clearly seen to couple effectively the 
singlet state with the £ + _ triplets through terms of the type J'S^S^ where I is 
one of the two sites of the dimer and m is a site of a neighboring dimer (actually 
the J' term couples coherently the triplets and the singlet on both neighboring 
dimers). When the singlet and the t + triplet are brought close to degeneracy by the 
field, this coupling becomes resonant and it changes drastically the property of the 
ground state. A critical field ~ J — aJ' (where a is a mo del- dependent constant) 
suffices to match the resonance conditions. The resonant coupling between the two 
states causes the breaking of all the symmetries previously retained by the reduced 
density matrix. A transfer of population from the singlet to the t + triplet breaks the 
symmetry between t + and £_, thereby leading to the appearence of a finite magnetization 
along the field. Moreover the reduced density matrix acquires finite off-diagonal terms 
(s|p( dim )|t+) = (t+|p( dim )|s) 7^ 0, which implies that, defining creation/annihilation 
operators for the t + bosons b\ \s) = \t+), b t+ \t + ) = \s), the expectation values (&] ) 
and (bt + ) become non-zero. This clearly identifies the breaking of the U(l) symmetry in 
the plane (i.e. the appearence of antiferromagnetic order transverse to the field) with 
the appearence of a condensate of t + bosons. 

From the magnetic point of view, the field-induced ordered ground state is a canted 
antiferromagnetic state of the type described at the beginning of the section, with finite 
uniform magnetization along the field and finite transverse staggered magnetization. 
By further increasing the field, the system is eventually brought to a fully polarized 
ferromagnetic state when reaching a critical field h& ~ J + aJ' (with a again model- 
dependent) which destroys the coherent mixing between the singlet and the t + triplet. 
Around this field the dimer reduced density matrix has p(\t + )) p(\s)), so the ordered 
state is conveniently described as a superfluid state of bosonic singlet holes in the 
triplet " sea" . The fully polarized state corresponds simply to the vacuum of such holes. 
Therefore we can conclude that a weakly coupled dimer system in a field realizes two 
successive insulator /superfluid transitions. 




Figure 1. (a) Heisenberg antifcrromagnet on a bilayer and in a magnetic field; a 
represents the layer index, (b) Same as in (a) but with site dilution. The open dots 
correspond to missing sites, and dashed J bonds correspond to dimers that have lost 
one spin. 
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3. Bilayer Heisenberg antiferromagnet in a field 

In what follows we consider the specific case of a S — 1/2 Heisenberg bilayer in a uniform 
magnetic field [Fig. Ufa)]: 

7~L = J' ' ~^ ^ ^2 ' ^ i ' 2 ~~ ^ €i ' a ^ta ■ (2) 

(ij) a=l,2 i i,a 

where the index i runs over the sites of a square lattice, (ij) are pairs of nearest 
neighbors on the square lattice, and a is the layer index. This model has been intensively 
investigated theoretically and it also captures the main features of the magnetic 

interactions in the recently synthetized compound BaCuSi206 [Ej. 

From the qualitative discussion of the previous section we have seen that the 
insulator/superfluid transition in a weakly-coupled dimer system involves essentially the 
singlet and one of the triplet states (t + ) of each dimer. Within a low-energy effective 
description, appropriate close to the degeneracy point between these two states, it is 
then reasonable to discard the other two triplet states (£+) and derive an effective 
dynamics for the triplet t + bosons only (or, equivalently, for the singlet bosonic holes). 
To retain the ground state symmetries of the dimer-singlet phase after discarding two 
of the triplets, it is necessary to assume that the ground state below the lower critical 
field h'S is a product of singlets \^f) = Yli \ s i), so that the bosonic insulating state is 
simply the vacuum for the t + triplets. 

The derivation of the effective bosonic Hamiltonian is well known ^H] and here 
we simply sketch it. Introducing the operators Sj = S^i + and tj = S^i — S it2 
the intradimer Hamiltonian reads 7imtra = (J/2) X]j( s i — ^ s f) (neglecting an irrelevant 
additive constant) and the interdimer one reads "Winter = (J'/fy J2{ij)( s i ' s j + U ' 
Truncating the dimer spectrum to \t + ) and \s), one introduces the pseudospin states 
t) = \t+), | I) = \s), and associated 5 = 1/2 pseudospin operators of, af. The 
Si, ti operators are then expressed in terms of the pseudospin ones as: sf = erf + 1/2, 
tf = \/2af, and all other components are zero within the subspace spanned by \s) and 
\t+). By further exactly mapping the pseudospins onto hardcore bosons, erf = — 1/2 
(n = tfb), af = b\, and = b i: with [bi, fej] = (i ^ j) and b\} = 0, we obtain the 
effective hard-core bosonic Hamiltonian 

Tucson = -j'/2 J2( b S + h - c -) + J, /2 J2 nin i - J ( h - x ) ni ( 3 ) 

(ij) {ij} i 

On each dimer location i, |0) = \s) and |1) = \t+). For h < hfi = J — zJ'/2 
the ground state is the vacuum, while for h > = J + zJ'/2 the ground state 
has one hardcore triplet per dimer (z = 4 is the coordination number of the square 
lattice). Both insulating states have a finite particle gap. For any intermediate 
field value < h < h!^ the system is gapless with bosons having a fractional 
filling (which corresponds to a finite uniform magnetization away from its saturation 
value, < m z u = (S*) < 1/2), and forming a condensate with long-range phase 
coherence (which corresponds to staggered antiferromagnetic order transverse to the 
field m x s {y) = ((-iysf y) )). 
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Fig. |21 shows the succession of ground-state bosonic phases in a bilayer Heisenberg 
antiferromagnet with g = J/ J' = 4. Shown are the results of a quantum Monte Carlo 
calculation based on the Stochastic Series Expansion method with directed-loop update 
[To] . The superfluid density T is estimated through winding number fluctuations [17], 
while the uniform susceptibility \u is obtained by numerical derivation of the uniform 
magnetization with respect to the field. 




0.5 1 1.5 2 2.5 



Figure 2. Field dependence of the uniform magnetization m*, uniform susceptibility 
Xu-, transverse staggered magnetization m x s , and superfluid density T at T = for the 
bilayer antiferromagnet with J / J' = 4 in the clean limit. 



4. Bilayer Heisenberg antiferromagnet with site dilution: Bose-glass phase 

A well controlled way of introducing lattice disorder in quantum spin systems is by 
doping the magnetic ions with non-magnetic ones with the same valence, so that the 
main effect is the removal of some of the spins in the magnetic Hamiltonian. In the case 
of BaCuSi 2 6 , e.g., this can be achieved by doping the Cu 2+ ions with non- magnetic 
ions as Zn 2+ or Mg 2+ . The magnetic Hamiltonian of the Heisenberg bilayer in presence 
of site dilution [Fig. UJb)] reads: 

TC — J ^Je^ie^S^l- <Si,2 + ^'^^ e i,a e j,a>Si,a '^j,a~Jh e i,ocSi,g -(4) 

i (ij) a=l,2 i,a 

Here the variables ej )Q , are random numbers taking value with probability p 
(corresponding to the concentration of non-magnetic dopants) and 1 with probability 
1 — p. 

In presence of site dilution, some of the dimer are completely missing, and some 
are reduced to single spins. We can nonetheless repeat the same approximate bosonic 
mapping as discussed in the previous section by truncating the Hilbert space of the intact 
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dimers while retaining the full space of the dangling spins. The J' coupling between two 
intact dimers has obviously the same expression in terms of bosonic operators as before. 
The coupling between a dimer and a dangling spin on the a-th layer is instead of the 
form Hdimer-spin = (J' /%) (*i + s i) ■ Sj jCt . Introducing the pseudospin operators for the 
dimers, we obtain ft dimer _ spin = (J'/y/i)(ofS? a + <T?Sl )0l ) + (J'/2)(<7? + 1/2)%. Finally, 
mapping the pseudospins onto b, hardcore bosons and the dangling 5 = 1/2 spins onto 
c, c* hardcore bosons, we end up with the following (rather complex) disordered bosonic 
model 



H 



boson 



Hhh + 7~t 



be 



(5) 



where 



n 



bb 



J' J' 



(ij) 



(ij) 



J 



1fi,ayj,c 



Mi Mi 



(ij),a 



n 



j' 



i+d 



+d) 



(ij) 



li Mi 



(7) 



be 



+ J' ^(AiTj ^M,- + 7jAj MiTij). 

(ij) 



(8) 



Here M- = eje 



T»,a = e»,a(l - £i,a), 7i = E a Ti,a, and A; = e^e^ and £ d runs over the 
four lattice vectors. 

For site dilution well below the bilayer percolation threshold (p < p* = 0.5244(2) 
[IB]), the above model contains a network of 6-sites (Aj = 1, 7« = 0), corresponding to 
intact dimers in the original spin Hamiltonian, interrupted by c-sites (A, = 0, ji — 1), 
corresponding to dangling spins in the original model, and by empty sites (A, = 0, 
7i = 0) corresponding to missing dimers. The most relevant feature of this complex 
model is that the chemical potential for the 6-bosons, fib = J(h — 1) is quite different 
from that of the c-bosons, \i c — Jh + y ^2 d (li+d — \+d)- In fact, if J 3> J 7 , the 
chemical potential for the c-bosons is always positive, and when hJ ~ J' the c-sites 
are almost all filled, which corresponds to field polarization of the dangling spins. On 
the contrary the chemical potential for the 6-bosons can be negative for h < 1, which 
means that, for the same field h, a system with only 6-sites can be empty of bosons, 
while a system with only c-sites would be completely filled. When b- and c-sites are 
combined together, c-bosons can hop to neighboring 6-sites, so that there is always a 
finite, albeit small, population of 6-bosons around c-sites, but, for negative enough 
bosons are prevented from traveling deep inside continuous regions of 6-sites. Therefore 
for h < hfi the bosonic model has all c-sites essentially filled and 6-sites only partially 
filled when neighboring a c-site, but continuous regions of 6-sites are essentially empty 
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of particles. If J >> J', in this situation the system has a gap to the addition of further 
particles, and it is in a quantum disordered state very similar to the insulating states in 
the clean limit. 

When h > hfj the first 6-bosons appear in a homogeneous system made exclusively 
of 6-sites, and they form a coherent condensate. In presence of disorder, instead, b- 
bosons are induced away from c-sites only in those rare regions which are locally well 
approximating a clean system, namely where the local coordination number is very close 
to z = 4 of the clean system. In regions with a lower local coordination number the 
local critical field is effectively higher, given its dependence on z. This simply means 
that the first bosons to appear in the bulk of the b regions are localized, although the 
6-sites might form a percolating cluster for low enough dilution of the lattice. Therefore 
we are in presence of a phenomenon of quantum localization of the 6-bosons, introducing 
a new phase in the disordered system, namely a Bose-glass phase. In such a phase there 
is no gap to the addition of a particle because, for any value of h > , we will find 
a region in the system where the local critical field equals h, so we can inject a boson 
there. This is therefore an unconventional quantum- disordered insulating phase with 
short-range correlations and a gapless spectrum. To overcome localization and enter 
the condensate phase, it is necessary to reach a higher critical field h c \ > hf^ at which 
a percolating network of 6-sites becomes accessible to bosons, so that long-range phase 
coherence can be established throughout the system. 

A similar phenomenological description can be repeated when approaching the 
upper critical field, h ~ . Here the description is simplified by the fact that, for 
such a high field, the c-sites can be thought of as perfectly filled, and bosonic holes are 
all on 6-sites (a bosonic hole on a c-site is energetically quite unfavorable). Therefore 
the c-sites can be simply quenched in the |Mj = 1) state, and we get the simplified 
bosonic model 



which is the hardcore Bose-Hubbard model on a site-diluted square lattice with random 
on-site chemical potential (fluctuations in the chemical potential appear only close to a 
site vacancy). As before, the regions of 6-sites with a lower local coordination number 
will have a local upper critical field which is less than hfj, which means that they 
tend to expel bosonic holes. For h < such holes remain localized on regions with 
local coordination close to that of the clean system, implying that, for a critical field 
h C 2 < hf-}-, the system loses the hole condensate and it enters a Bose-glass phase of 
bosonic holes. 

5. Diluted Heisenberg bilayer: quantum Monte Carlo data 



7~Lb< 



■oson 




\ J{h-l)-J'Y,H+d (9) 



i L d 



The approximate bosonic picture we put forward in the previous section for the physics 
of the site-diluted Heisenberg bilayer is confirmed by extensive quantum Monte Carlo 
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simulations done on the original spin model making. Making use of Stochastic-Series- 
Expansion quantum Monte Carlo, we were able to efficiently reach the physical T = 
behavior by a successive increase of the inverse temperature (THj. We investigated system 
sizes up to 40 x 40 x 2, averaging typically over 200 disorder realizations. 



Figure 3. Zero-temperature field scan in the site-diluted bilayer Heisenberg 
antifcrromagnet with J/ J' = 4 and p = 0.2. Inset: scaling of the order parameter 
m x a in the different phases of the system; from top to bottom h = 0.75 (XY ordered), 
0.05 (order-by-disorder), 1.85 (high-field Bose glass), 0.35 (quantum disordered), 0.6 
(low-field Bose glass). 



Fig. H2 shows the succession of phases in a bilayer Heisenberg model with J / J' = 4 
and site dilution p = 0.2 when the applied field h is changed from zero to the saturation 
value h = 2. This picture has to be contrasted with the analogous field scan in the clean 
limit, as shown in Fig. 121 

At low field, a particular feature of the spin model, not discussed previously in 
the bosonic mapping, shows up. Doping the dimer-singlet phase, we get long-range 
antiferromagnetic ordering of the free moments appearing around each of the vacancies, 
a well-known phenomenon of order-by-disorder This corresponds approximately 

to bosons being created around c-sites but largely fluctuating in number given that the 
chemical potential \l c is very small, and coherently propagating between c-sites through 
effective long-range hoppings which decay exponentially with the inter-site distance, 
so that tenuous long-range phase coherence is established. Increasing the field leads 
to a strong increase in the fi c chemical potential, which quenches particle number 
fluctuations and it destroys the disorder- induced superfluid order. The system enters 
then a gapless quantum-disordered phase, in which the gapless nature, reflected in the 
finite compressibility of the bosons (susceptibility of the magnetic system), is due to 
the fact that c-sites are not all completely filled. For a larger ratio J / J' than the one 
considered here we can access the regime where c-sites are essentially all filled and the 
system does not admit any further particle, namely the compressibility vanishes and the 




\ • L=24 
• ■ L=32 
\ ♦ L=40 



h 
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system acquires a gap as in the clean case [21]. In this regime, either gapful or gapless, 
bosons are mostly pinned to c-sites, and they cannot propagate in 6-site regions due to 
the local finite particle gap in such regions. 




• l-m = 0.1 • l-m = 0.1 

(b) 

Figure 4. (a) Real-space images of the dimer magnetization rrij = (Sf j + Sf 2 ) on 
intact dimers in a 40x40x2 bilayer with J j J' — 4, dilution p = 0.1 and at inverse 
temperature j3J = 256, for h = 0.56 (left) and h = 0.6 (right). The radius of the 
dots is proportional to the dimer magnetization. The magnetization of unpaired spins 
(c-sites) is omitted for clarity. The most visible localized states are highlighted in the 
left panel, while the backbone of the percolating magnetized network is highlighted 
in the right one. The central panel shows the superfluid density as a function of the 
field for the specific sample considered, (b) Real-space images of the distance of the 
dimer magnetization from its saturation value, 1 — m^, for h — 1.87 (left) and h — 1.93 
(right). Other symbols, parameters and explanations as in (a). 

When h > ~ 0.47(1), nonetheless, large clean 6-regions accept the appearence 
of bosons, which remain quantum-localized, giving rise to a Bose-glass phase. This 
is clearly seen in real-space picture of the dimer magnetization m ; = (S-^ + S? 2 ) in 
Fig. Efa) (left panel), representing the local density of bosons on 6-sites within the 
approximate bosonic mapping. We observe that, in this regime, the order parameter 
m x s and the superfluid density are zero, and the number of bosons, corresponding to the 
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uniform magnetization, increases extremely slowly with increasing field, corresponding 
to the gradual appearence of particles in the rare clean regions of the system. 

Increasing the field beyond h c i ~ 0.69(2) a superfluid state is finally established 
in the system through derealization of the collective state of the bosons [Fig. Efa) 
(right panel)]. Switching then to a description in terms of bosonic holes when the 
filling exceeds 1/2, we observe that the delocalized superfluid state of holes, whose 
local concentration is proportional to 1 — mi [Fig. Efb) (left panel)], is destroyed at 
a critical field h C 2 ~ 1.78(2) lower than the saturation one through a mechanism of 
quantum localization [Fig. Hfb) (right panel)], which gives rise to the hole Bose-glass 
phase at high fields. Even more clearly than in the lower-field case, this phase shows 
a finite compressibility and hence a gapless spectrum, but absence of long-range phase 
coherence due to quantum localization. 

6. Conclusions 

In this work we have investigated the emergence of a rich physical scenario of 
strongly correlated bosons in weakly-coupled dimer antiferromagnets in a magnetic 
field. In presence of site dilution of the magnetic lattice, extended Bose-glass 
phases appear around the field-induced ordered phase, suggesting the possibility of an 
experimental realization of such a phase in quantum magnets. The particular case 
of a bilayer antiferromagnet we considered in this paper is realized by BaCuSi2C>6 
|14j . Nonetheless the picture of quantum localization of bosonic quasiparticles 
applies straighforwardly to other spin gap systems with different geometries, such as 
Sr 2 Cu(B0 3 )2 [22| or Tl(K)CuCl 3 jl], where site dilution can be analogously realized 
by doping the Cu sites. The availability of high magnetic fields allows to scan the 
whole succession of phases discussed in the paper, and different probes, such as elastic 
neutron scattering and magnetometry measurements, can directly access the magnetic 
observables characterizing the various phases. An open question is how to directly probe 
the superfluid nature of the bosonic quasiparticles condensed in the ordered-ground 
state, given that a supercurrent of such particles corresponds to a pure spin current 
with no energy or charge transport associated with it. 
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